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Abstract 

We study the dynamical transport in weakly coupled superlat- 
tices in the presence of intense radiation in the THz regime. We 
derive a general model for the time dependent tunneling current 
within the Keldysh nonequilibrium-Green-function formalism. 
For the particular situation in which fast scattering procceses 
drive the system to local equilibrium in the wells drastic changes 
are found in the current vs. voltage curves. 

Interaction with external time-dependent fields in semiconductor superlattices (SLs) leads 
in many cases to completely new ways of electrical transport: negative pumping of electrons, 
dynamical localization, photoinduced electric field domains (EFDs) , and bistability between 
positive and negative current [1,2]. In voltage biased weakly-coupled SLs, stationary EFDs 
arise if the well doping is large enough. They usually consist of two homogenous electric 
field regions separated by a domain wall (DW) of accumulated electrons. When the carrier 
density decreases below a critical value, the DW is unable to find a stable position and 
starts to move over several periods. Its motion and recycling gives rise to spontaneous self- 
oscillations of the electric current [3]. As doping is experimentally hard to manipulate, this 
static-to-dynamic-state transition has been achieved by the aplication of transverse magnetic 
fields [4], laser illumination in undoped SLs [5], and careful variation of temperature [6]. 

Due to the strongly nonlinear behaviour of a weakly-coupled SL (stemming from the 
interplay between tunneling processes and averaged Coulomb interactions), perturbing the 
system with an ac potential involves the rising of even more complex aspects in the physics 
of transport. In particular, it has been reported analyses of time-dependent current in the 
presence of a low-frequency signal both experimentally and theoretically [7]. In that case, 
the ac frequency is of the order of tens of MHz and the ac potential results in an adiabatic 
modulation of the system. In this paper we are interested in a very different regime. We 
investigate the time-dependent current through a multiple quantum well driven by a high- 
frequency ac potential, V°~ c {t) = V" 10 cos(cuot) , where the ac frequency f ac = oj^jlix is of the 
order of several THz and V^ c is the ac amplitude in the zth quantum well (QW). It is well 
known that in this case photoassisted tunneling takes place, and the electronic states develop 
side-bands which act as new tunneling channels. 

The Keldysh Green-function formalism [8] allows to obtain general expressions for the 
tunneling current. It also permits to include electron-electron interaction in the system. 
Both the limits of noninteracting systems and local equilibrium can be deduced within this 
scheme. We shall restrict ourselves to analyze the case of a weakly-coupled SL where the 
electron-electron interaction will be included in a mean-field (Hartree) manner. From the 
time evolution of the occupation number operator of the ith QW (iVj), the change in the 
number of electrons in the ith QW is e(Ni) = /j-i,, — This continuity equation relates 

the change of the number of electrons in the ith QW with the current density flowing from 



the i — 1st (ith) QW to the ith (i + 1st) QW, Ij-i^ (J^+i). By using the standard techniques 
of motion equation [8] one arrives to 



Am (*) = | Re E T Ui+i J dr [Gl +i (t, T )gl (r, t) + (t, r)<?| (r, t) 
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x being the transmission coefficient, ki represents the set of quantum numbers which 

label the electronic states in the ith QW. Here, g^"^ is the advanced (lesser) Green function 

(GF) for a decoupled QW in the presence of an ac potential and scattering processes. GC 
is the retarded (lesser) QW GF when scattering, ac signal and tunneling are present. In the 
case of a SL, the GFs which appear in Eq. ([!]) can be written down exactly by making several 
assumptions. Firstly, one considers the scattering self-energy (due to, e.g., impurities or LO 
phonons) in a phenomenological way as an energy indenpendent constant (which is denoted 
here by 7 = Im£ sc ). As in weakly-coupled SL's the scattering lifetime (~ 1 ps) is much 
shorter than the tunneling time (~ 1 ns), we can assume local equilibrium within each QW 
and neglect the tunneling self-energy in the expressions of the GFs. The electron transport is 
thus sequential. Once these approximations are made, the expression for the retarded (lesser) 

ieV ac 

ac-driven QW GF in the presence of scattering is G^ (<) (t, t') = e^^^-^^Q^ (t-f), 

where G\ (t — t') is the static retarded (lesser) QW GF's in the presence of scattering. 
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local equilibrium reads [9] 

2e m=oo , 

Ji,i+i{t) = -z- T k t k l+1 E Jm(Pi,i+i)\ cos ([3 iji+ i sin u t -mw t) (2) 
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Similar relations hold for gi '(t,t'). By using these expressions the current in the case of 
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A k i+1 ( e + mhuj )A^(e) (fi(e) - f i+1 (e + mhu ))j + sin(/3 iji+ i sin (u t - muj t) 
A k 1+1 ( e + m ^o)Re^|^(e)/ i+ i(e + mtkv ) + Re(5^. +i g. +i (e + mfuv )A^(e)fi(e) 

where v4g denotes the spectral function for the ith isolated QW in the presence of scattering. 
J p is the pih Bessel function whose argument is given by (3 = e(Vf c — V^^/fiujQ and fi(e) is 
the Fermi-Dirac distribution function for the ith QW. 

Notice that the current may be written as I(t) = Io + J2i>o If° s cos (^0 t) +If m sin {lui t), 
where Iq is the time-averaged current [1]. I™ 3 and If m contain higher harmonics for / > 
0. Now, since we are interested in the photoassisted current, ujq is much larger than the 
tunneling rate. In addition, the scattering lifetime represents the lowest temporal cutoff 
above which our assumption of local equilibrium within each QW holds. In other words, 
to ensure the vality of Eq. (|3|) one is restricted to study dynamical processes whose time 
variation is longer than h/j. As intense THz fields typically have hu > 7 we must carry 
out a time average of I(t). This implies that the explicit time variation of I(t) vanishes and 
we are left with the implicit change of Iq with respect to time. This variation (in scales 
larger than h/j) results from the evaluation of the continuity equation for i = 1,...,N, 
where N is the number of wells, supplemented with Poisson equations, constitutive relations, 
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Figure 1: (a) Time-averaged current vs. ap- Figure 2: (a) Time-averaged current vs. ap- 
plied voltage (for parameters see text), (b) plied voltage as (3 is varied, (b) I(t) for 
Time dependence of the current for V=1.1V V = 1.1 V. Curves for (3 ^ are shifted for 
for the values of (3 shown in (a). clarity. 



and appropiate boundary conditions (see [9] for details). Of course, by incorporating an 
accurate microscopic model for the scattering dynamics and calculating the nonequilibrium 
distribution function for each QW in the presence of the ac signal, one could investigate all 
the parameter range. This is out of the scope of the present work. 

The total current traversing the sample consists of the tunneling plus displacement cur- 
rents, I(t) = JiA + i + (e/d)(dVi/dt), where (e is the static permittivity, d the barrier width, 
and V% the voltage drop in the ith barrier. Given reliable initial conditions, If) is self- 
consist ently calculated for a N = 40 SL with 13.3-nm GaAs wells and 2.7-nm AlAs barriers. 
Well doping is 2 x 10~ 10 cm -2 and we take 7 = 8 meV and f ac = 3 THz. In Fig. |T](a) 
the average of If) is plotted as a function of the applied dc bias, V. Without ac forcing 
the I— V curve shows branches after the peak corresponding to resonant tunneling between 
the lowest subbands (see inset). This feature is distinctive of static EFD formation [1-3]. 
In the presence of an ac signal the branches coalesce and a plateau is formed- this is the 
key signature of current self-oscillations. I(t) for V = 1.1 V is depicted in Fig. |l|(b). For 
(3 = the current achieves a constant value after a transient time. As (3 is increased damped 
oscillations are observed until a stable current self-oscillations arise when (3 = 1. This is an 
indication that the ac potential induces a transition from a stationary configuration towards 
a dynamic state possibly via a supercritical Hopf bifurcation. The oscillation frequency is 
of the order of 150 MHz, much smaller than uq. Then we conclude that the existence of 
sidebands and its influence over the non-linear behaviour of the system drives the SL towards 
oscillations. 

Increasing 7 results in dramatic consequences. This may be effectively achieved by either 
applying a transverse magnetic field [4] or raising the temperature [6]. In the absence of ac 
potentials for 7 = 11 meV there is a voltage range (~0.8-1.3 V) where now self-oscillations 
show up [3] (see Fig. @(a)). As (3 enhances the plateau starts to be replaced by a positive 
differential resistance region. There is a similar well-known phenomenon in weakly-coupled 
SLs: under a critical value of the carrier density neither static nor moving DWs exist and the 
electric field drops homogenously across the whole sample. In our case the doping density 



is constant and it is the ac field the physical parameter which forces this transition. To 
illustrate this we have calculated I(t) for a fixed bias 7= 1.1 V (see Fig. 0(b)). At (3 = 
the current oscillates with a frequency of 170 MHz. This is a result of the motion of the 
accumulating layer of electrons, and its recycling in the highly-doped contacts [3]. For (3 = 2 
the current is already damped, and I(t) reaches rapidly a uniform value for (3 = 4. This 
is a striking feature- an oscilation disappearance induced by an ac potential. A qualitative 
explanation follows. 

The general transition (static EFDs)— ^(moving EFDs) -^(homogenous electric field pro- 
file) depends basically on the particular shape of the drift velocity between two wells [3-6]. 
For a sufficiently high peak-to- valley ratio the system evolves towards a static EFD situation. 
If this ratio is very low the applied voltage drops smoothly across the sample. In between, 
self-sustained current oscillations show up. In our case the application of a THz field pro- 
duces tunneling with absorption and emission and photons, thus strongly modifying the 
current between two wells (or, roughly, the drift velocity). In particular, for the values used 
in the text the zero-photon tunneling peak (weighted by Jq{0)) is depressed and the current 
associated with the emission of, e.g., one photon becomes comparable when (3 increases from 
0. Then the peak-to- valley ratio decreases as (3 increases and the current oscillations are 
reached provided we start from a stationary current (see Fig. [I]). On the other hand, if the 
dynamic configuration is stable for (3 = the ac potential will tend to drive the SL to a 
trivially homogenous electric field profile (see Fig. 0). A more detailed analysis will be given 
elsewhere [9]. 
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